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ABSTRACT 

Using an % — z dropout criterion, we determine the space density of z ~ 6 
galaxies from two deep ACS GTO fields with deep optical-IR imaging. A total 
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of 23 objects are found over 46 arcmin 2 , or ~0.5±0.1 objects arcmin -2 down to 
zab ~ 27.3 (6cr), or a completeness-corrected ~0.5±0.2 objects arcmin -2 down 
to zab ~ 26.5 (including one probable z ~ 6 AGN). Combining deep ISAAC 
data for our RDCS1252-2927 field (J AB ~ 25.7 and Ks AB ~ 25.0 (5 a)) and 
NICMOS data for the HDF North (J 110 ,ab and H wo ,ab ~ 27.3 (5 a)), we ver- 
ify that these dropouts have relatively flat spectral slopes as one would expect 
for star-forming objects at z ~ 6. Compared to the average-color {(3 = —1.3) 
[/-dropout in the Steidel et al. (1999) z ~ 3 sample, i-dropouts in our sample 
range in luminosity from ~1.5 (zab ~ 25.6) to ~0.3 L* (zab ~ 27.3) with 
the exception of one very bright candidate at Zsso.ab ~ 24.2. The half-light radii 
vary from 0.09" to 0.29", or 0.5 kpc to 1.7 kpc. We derive the z ~ 6 rest-frame 
UV luminosity density (or star formation rate density) using three different pro- 
cedures. All three procedures make use of simulations based on a slightly lower 
redshift [z ~ 5) V^-dropout sample from CDF South ACS images. First, we 
make a direct comparison of our findings with a no-evolution projection of this 
^-dropout sample, allowing us to automatically correct for the light lost at faint 
magnitudes or lower surface brightnesses. We find 23±25% more i-dropouts than 
we predict, consistent with no strong evolution over this redshift range. Adopting 
previous results to z ~ 5 (Bouwens, Broadhurst, & Illingworth 2003; Thompson 
et al. 2001), this works out to a mere 20 ± 29% drop in the luminosity density 
from z ~ 3 to z ~ 6. Second, we use the same ^/-dropout simulations to derive 
a detailed selection function for our z-dropout sample and compute the UV- 
luminosity density (7.2 ± 2.5 x 10 25 ergss -1 Hz -1 Mpc -3 down to zab ~ 27). We 
find a 39 ±21% drop over the same redshift range [z ~ 3 — 6), consistent with the 
first estimate. This is our preferred value and suggests a star formation rate of 
0.0090 ± 0.0031 M Q yr -1 Mpc" 3 to z AB ~ 27, or ~ 0.036 ± 0.012 M yr -1 Mpc" 3 
extrapolating the luminosity function to the faint limit. Third, we follow a very 
similar procedure, except that we assume no incompleteness, finding a rest-frame 
continuum luminosity density which is ~ 2 — 3 x lower than our other two deter- 
minations. This final estimate is to be taken as a lower limit, and is important 
in the event that there are modest changes in the colors or surface brightnesses 
from z ~ 5 to z ~ 6 (the other estimates assume no large changes in the in- 
trinsic selectability of objects). We note that all three estimates are well within 
the canonical range of luminosity densities (e.g., Madau, Haardt, & Rees 1999) 
necessary for reionization of the universe at this epoch by star-forming galaxies. 

Subject headings: galaxies: evolution — galaxies: formation - galaxies: high- 
redshift 
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1. Introduction 

The Hubble Deep Field (HDF) campaign has been highly influential in shaping our 
understanding of star formation in the high-redshift universe (Williams et al. 1996; Casertano 
et al. 2000; Ferguson, Dickinson, & Williams 2000). Early results demonstrated that the star 
formation rate density-as measured from the rest-frame continuum UV-increased from z ~ 4 
to an apparent peak around z ~ 1 — 3 (Lilly et al. 1996; Madau et al. 1996; Connolly et 
al. 1997; Cowie et al. 1999). While these results were largely solidified by Steidel et al. 
(1999) with his wide-area [/-dropout survey, the addition of NICMOS data to the HDF 
North demonstrated that this trend continued to z ~ 6 (Thompson et al. 1999; Deltorn et 
al. 2003; Bouwens, Broadhurst, & Illingworth 2003, hereafter denoted BBI). 

Unfortunately, these studies were limited enough in area to raise doubts about how 
representative they really were of the high-redshift universe. They also suffered from a lack 
of deep high-resolution imaging at wavelengths intermediate between the optical regime and 
the near infrared, necessary for obtaining a more detailed look at the z ~ 6 — 7 universe. 
Deep NICMOS observations have been useful in addressing this latter shortcoming, but 
only partially due to its small field of view. Fortunately, the installation of the Advanced 
Camera for Surveys (ACS) (Ford et al. 1998) on the Hubble Space Telescope has helped to 
redress several of these issues, including crucially for the first time imaging in the z-band, 
permitting a more secure detection of objects at high redshift (z ~ 5 — 6). Moreover, its 
10 x improvement over WFPC2 in surveying capability allows large areas to be surveyed 
to nearly HDF depths (Ford et al. 2003), the Great Observatories Origins Deep Survey 
(GOODS) being a notable example (Dickinson & Giavalisco 2002). 

Here, we describe some early work done using deep ACS data to extend these searches 
to z ~ 6, to establish the prevalence of galaxies in this era. Interest in star formation at 
z ~ 6 has been particularly intense as of late because of recent absorption line studies on 
3 QSOs at z > 5.8, suggesting that reionization may have happened at about this epoch 
(Fan et al. 2001; Becker et al. 2001; Fan et al. 2002). In this work, we consider two fields 
from the ACS GTO program, RDCS 1252-2927 and the HDF North, in our search for z ~ 6 
objects. Both fields have deep ACS i and z data, and infrared observations, important for 
securely identifying z ~ 6 objects. Relative to other work (Yan et al. 2003; Stanway, Bunker, 
Sz McMahon 2003), the present search is slightly deeper, with better IR data to confirm the 
redshift identifications. In fact, our use of the HDF North field is especially propitious, 
given the exceptionally deep WFPC2 and NICMOS images available to examine faint z ~ 6 
candidates. We put this new population in context by comparing them with lower-redshift 
expectations, projecting z ~ 5 galaxy samples from CDF South GOODS to z ~ 6 using 
our cloning formalism previously used in work on the HDFs (Bouwens, Broadhurst, & Silk 
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1998a,b; BBI). 

We begin by presenting our data sets, describing our procedure for doing object detection 
and photometry, and finally discussing our z ~ 6 z-dropout selection criterion (§2). In §3, we 
present our results. In §4, we describe a comparison against the wide-area GOODS sample 
and then use these simulations to make three different estimates of the z ~ 6 rest-frame 
continuum UV luminosity density (§5). Finally, in §6 and §7, we discuss and summarize our 
findings. Note that we denote the F775W, F850LP, F110W, and F160W bands as i 775 , 
z S50 , J 110 , and #160) respectively, and we assume [Qm^aM = [0.3,0.7,0.7] in accordance with 
the recent Wilkinson Microwave Anisotropy Probe (WMAP) results (Bennett et al. 2003). 

2. Observations 
2.1. Data 

There are two different fields from our ACS GTO program which are particularly useful 
for 2-dropout searches. The first involves deep ACS WFC i 775 and z 8 50 images of RDCS1252- 
2927, a z = 1.235 cluster. RDCS1252-2927 was selected from the ROSAT Deep Cluster 
Survey (Rosati et al. 1998; Rosati et al. 2003). Three orbits in i 775 and five orbits in z 850 
were obtained at four overlapping pointings, arranged in a 2x2 grid with an overlap of ~1 
arcmin so that the overlapping regions (~10 arcmin 2 ) were covered to a depth of six orbits 
in i 775 /ten orbits in z 850 with a small central region (~1 arcmin 2 ) being covered to a depth of 
twelve orbits in 2 775 /twenty orbits in zgso- The ACS images were aligned, cosmic-ray rejected, 
and drizzled together using the ACS GTO pipeline (Blakeslee et al. 2003a). 

Very deep integrations were obtained on ISAAC over 4 overlapping regions on RDCS1252- 
2927 (covering 4x4 arcmin, or ~44% of our 36 arcmin 2 ACS mosaic). A total of 24.1 and 22.7 
hours were invested in the J and Ks integrations, respectively (~6 hrs and ~5.8 hrs at each 
of the 4 offset positions). These observations reached Jab = 25.7 (5a) and Ksab = 25.0 
(5a) in the shallower, non-overlapping regions and Jab = 26.5 (5<r) and Ksab = 25.8 (5a) 
in the small (1' x 1') central region, with a FWHM for the PSF which was almost uniformly 
~ 0.45" across the entire IR mosaic. These data were then aligned with our optical data 
and resampled onto the same 0.05"-pixel grid. 

The second field utilizes deep ACS observations of the HDF North, taken as part of our 
GTO program, 2.5 orbits in i 775 and 4.5 orbits in z S50 . This data is supplemented with 1.5 
orbits in i 775 and three orbits in 2 850 from the GOODS program in this field (representing 
three epochs of the GOODS program) to yield a total depth of 4 orbits in ^775 and 7.5 orbits 
in ^850 over an effective area of 10 arcmin 2 . To complement the ACS % and z data, both the 
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HDF North optical data (Williams et al. 1996) and JH infrared data from the Dickinson 
(1999) campaign were aligned and registered onto the same 0.05"-pixel scale as our ACS 
fields, leaving the WFPC2 data with a FWHM of ~ 0.18" for the PSF and the NICMOS 
data with a FWHM of ~ 0.25" . The NICMOS images reached J 110 ,ab ~ 27.3 (5<r) and 
H 160 ,ab ~ 27.3 (5a). 

The i 775 ,AB = 25.64 and z 850AB = 24.84 CALACS (02/20/03) zeropoints (Siranni et 
al. 2003) are assumed throughout, along with a galactic absorption of E(B — V)=0.075 
and 0.012 for the two fields (from the Schlegel, Finkbeiner, & Davis 1998 extinction maps), 
resulting in a correction of -0.11 m and -0.02 m to the zsso zeropoint for RDCS1252-2927 and 
the HDF North, respectively (and a correction of -0.15 m and -0.024 m for the ^775 filter). 

2.2. Detection and Photometry 

Briefly, object detection is performed on the basis of our deep WFC z 850 images after 
smoothing the images with a 0.09"-FWHM Gaussian kernel and looking for 4.5<r peaks. 
Photometry is obtained for all detected objects with SExtractor (Bertin & Arnouts 1996) 
using two scalable apertures, the inner one to measure colors and the outer aperture to 
estimate the total flux. For both sets of IR data, similarly-scaled apertures are used to 
measure colors, but a correction is applied based on the z 850 image to estimate how much 
flux is lost in the IR due to PSF smoothing. Due to correlation in the noise, a concerted 
attempt was made to model the noise so that a relatively realistic treatment of uncertainties 
could be applied throughout the analysis (Appendix A). A more comprehensive description 
of our techniques for object selection and photometry is given in BBI. 16 

In total, 4632 and 1261 objects were recovered in the RDCS1252-2927 and HDF North 
fields, respectively. Most spurious detections were eliminated by demanding that each object 
be a 6<T detection within one Kron-radius (Kron 1980) (typically ~ 0.15"). Areas contam- 
inated by optical ghosts or satellite trails were not included in the analysis (the excluded 
area was < 0.5%). A number of spurious detections were found around bright stars or el- 
lipticals, a problem exacerbated by the rather extended wings on the z 8 5o PSF. After some 
preliminary cleaning, all point-sources were removed from our catalogs (~ 503 and ~ 110), 
the rationale being to eliminate very red stars which might otherwise masquerade as high- 
redshift objects. We found that the SExtractor stellarity parameter adequately identified 
stellar objects. While such a cut might eliminate genuine z ~ 6 star-forming objects, all 



16 Notc that our procedure for object detection and photometry is different from that used by the GTO 
team (Blakeslee et al. 2003a; Benitez et al. 2003). 
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of the very red ((2775 — z 8 ^o)ab > 1.5) point-like objects (~ 4) were found to have much 
redder z — J colors (~ 1.2) than most of the probable z ~ 6 objects, and therefore a stellar 
identification seemed reasonable (see the systematic color differences between point-like and 
extended (i 775 — z 850 )ab > 1-5 objects presented in Figure 1). We note that an examination 
of red i — z > 1.3 point-like objects in our fields revealed one probable z ~ 6 AGN (Appendix 
B). 



2.3. Dropout Selection 

The Lyman-break technique takes advantage of the increasingly strong deficit of flux at 
high redshift caused by the intervening Lyman-alpha forest eating into the spectrum short- 
ward of 1216 A (Madau 1995). Combining this with flux information redward of the break 
permits one to determine the spectral slope redward of the break and therefore relatively 
robustly distinguish the objects of interest from intrinsically red galaxies at lower redshift, as 
demonstrated by extensive spectroscopic work done on a variety of different dropout samples 
(Steidel et al. 1996a,b; Steidel et al. 1999; Weymann et al. 1998; Fan et al. 2001). 

For our filter set, the i 77 s — z$5o color measures the spectral break and the z 8 ^q — J 
color defines the spectral slope redward of this break. In Figures 1-2, we illustrate how a 
starburst spectrum (100 Myr continuous star formation) attenuated with various opacities 
of dust (E(B — V) = 0, 0.2, 0.4) would move through this color-color space as a function of 
redshift. In both plots, it is clear that beyond z ~ 5.5, the template i — z colors become 
very red (> 1.2) while the z — J colors remain very blue (< 0.5 — 1). For reference, we also 
include the colors of possible lower redshift interlopers using the Coleman, Wu, & Weedman 
(1976) spectral templates. 

After considerable experimentation, we adopted a simple (2775 — z 8 5o)ab > 1-5 cut 
throughout in our selection of z-dropouts. This color cut is motivated by evaluating object 
selection in regions where we have both ACS and infrared data. We present such data 
in Figures 1-2 from our RDCS1252-2927 and HDF North fields, providing 2a upper-limits 
for objects without significant infrared flux. We have also lightly shaded those regions in 
i — z, z — J color space where z > 5.5 i-dropouts are expected to lie: (2775 — z$5o)ab > 1-5, 

- Zb6o)ab > (z 775 - J) AB + 0.7, (^ 775 - J)ab < 0.9, z 850AB < 27.3 for RDCS1252-2927 
and (i 775 - z 850 ) AB > 1.5, (i 775 - z 850 )ab > (z 775 - J 1w )ab + 1-0, (2775 - J\xo)ab < 0.6, 
Z85q,ab < 27.3 for the HDF North. A quick glance shows that objects with very red (> 1.5) 
i — z colors also have blue z — J colors and lie exclusively in this region, thereby validating 
our basic selection criteria. Over the 21 arcmin 2 where we have infrared coverage, we find 11 
objects in RDCS1252-2927 and 1 object in the HDF North which satisfy our i — z > 1.5 cut. 
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Fig. 1. — (z 775 — z$5o)ab vs. (-2850 — J) ab colour-colour diagram illustrating the position of 
our RDCS1252-2927 i-dropout sample (shaded region) relative to the photometric sample 
as a whole. Tracks for a 10 8 year starburst with various amounts of extinctions have been 
included to illustrate both the typical redshifts (labelled for z = 5.5 and z = 6) and SED 
types included in the selection window. The low-redshift (0 < z < 1.2) tracks for typical E, 
Sbc, and Irr spectra have been included as well to illustrate the region in colour-colour space 
where possible contaminants might lie. There is a clear separation between the % — z > 1.3 
point-like objects (crosses) and i — z > 1.3 extended objects (solid squares) along the z — J 
axis. The distribution of objects in color-color space led us to adopt (i — z) > 1.5 as our 
generalized i-dropout selection criteria. In all cases, error bars represent 2 a limits. The 
clump at (z 775 — z 850 )ab ~ 0.9 and (^ 850 — J)ab ~ 1 ar e early-type galaxies from the cluster. 
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Fig. 2. — Same as Figure 1, but for i-dropouts from the HDF North. Note that the NICMOS 
Jiw,ab filter shown here is distinct from the ground-based Jab filter used in Figure 1. 
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^775^850 J K s and ^775^850^110-^160 images for these objects are shown in Figure 3, along with 
plots showing their position in color-color space, fits of plausible SEDs to the broadband 
fluxes, and an estimated redshift. The photometric redshifts are estimated using a bayesian 
formalism similar to that outlined in Benitez (2000), using a prior which matches the observed 
distribution of z ~ 3 spectral slopes (Steidel et al. 1999). Our typical 2-dropout has a S/N 
of ~ 2 — 3 in the infrared, more than adequate to put good constraints on the spectral slope 
redward of the break. To illustrate this and to indicate how different the i-dropouts really are 
from possible low redshift contaminants, we show three i — z > 1.3 objects with red z — J > 0.8 
colors in Figure 4. For reference, we also include a figure with the z = 5.60 Weymann et 
al. (1998) object from the HDF North to show the position of a spectroscopically-confirmed 
z > 5.5 object on these diagrams (Figure 5). 

Having presented our entire sample of z-dropouts with infrared coverage, we now move 
onto quantifying the contamination rate due to low redshift interlopers. The most obvious 
way of doing this is simply to count the fraction of objects with i — z > 1.5 which satisfy 
the two-color criteria we specified above versus those that do not. Unfortunately, for many 
objects we only have limits and not precise measures of the IR colors, leaving us with cases 
where we are not sure if an object lies in our sample or not. Therefore, following Pozzetti et al. 
(1998) in their analysis of dropouts in the HDF North, we resort to the use of Kaplan-Meier 
estimators with censoring (Lavalley, Isobe, & Feigelson 1992), where the implicit assumption 
is that censoring is random, e.g., that objects with and without limits are drawn from the 
same parent distribution. Given the narrow range of magnitudes and sizes in our sample we 
believe this assumption to be approximately satisfied. Performing this analysis on all objects 
in our fields with % — z > 1.3 (46 objects), we find the following contamination fraction as 
a function of (i - z) color: 0% (i - z) > 1.7, 13% (i - z) > 1.5, and 21% («-*)> 1.3. 
(Note that this estimate is for samples from which the point sources (stars or AGNs) have 
already been removed.) While our data set contains two % — z ~ 1.3 — 1.4 objects which are 
apparently ellipticals or Extremely Red Objects (EROs) at ~25.3, z — K ~ 2 and z — K ~ 3 
(Figure 4), all the objects with (i — z) > 1.5 have infrared colors consistent with their being 
at high redshift. 

This small contamination fraction (13%) allows us to substantially increase the size of 
our z ~ 6 sample by including (z 775 — z 850 )ab > 1-5 objects without IR coverage. A total of 
11 such objects satisfy these criteria: 10 from RDCS 1252-2927 and 1 from the HDF North. 
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Fig. 3a. — in5Z 850 JK postage stamps of the 12 z-dropouts (z ~ 6) identified in our ACS 
fields over the 21 arcmin 2 where we had infrared coverage (11 z-dropouts did not have 
measurements in the IR). ^775^850^110-^160 postage stamps are shown for the object from 
the HDF North. The optical and IR images are smoothed with 3x3 and 6x6 boxcars, 
respectively. The position in % — z, z — J space is also indicated along with the broadband 
SEDs and estimated redshift. As in Figures 1-2, we have included lines denoting the way 
starburst objects (10 8 yr bursts) with various dust attentuations would move through color- 
color space as a function of redshift. We have also included the tracks of possible interlopers. 
As in Table 1, the "1252-" prefix denotes an object from RDCS1252-2927. The postage 
stamps are 3.0" x 3.0" in size. 
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Fig. 3b.- 




Fig. 3c— 
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Fig. 4 

• — ^775^850^-^ postage stamps of several i — z > 1.3 objects with moderately red 
z — J colors. In all three cases, objects are clearly visible in the infrared, illustrating the 
basic value of our infrared data for distinguishing high-redshift objects from low redshift 
interlopers. The top two objects are likely to be low- redshift (z ~ 1.2 — 1.5) ellipticals. The 
third object, though likely at high redshift, did not make our z-dropout selection cut. The 
postage stamps are 3.0" x 3.0" in size. 
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Fig. 5. — ^775^850<^iio-^no postage stamps of the z = 5.60 Weymann et al. (1998) object 
from the HDF North (HDF4-473.0). While this object did not meet our selection criteria, 
it is reassuring to find a spectroscopically-confirmed object at the low-end of our i-dropout 
selection range (see §4) just missing our (i 775 — z S50 )ab > 1-5 cut to the blue. We measured 
an (z 775 — z 850 )ab = 1-2 color for this object. 



3. Results 



In summary, we find 21 objects over 36 arcmin 2 in our RDCS1252-2927 field and 2 
objects over 10 arcmin 2 in our HDF North field which satisfy the z-dropout criterion (i — z > 
1.5) we defined in the previous section. There is deep IR coverage for 12 of these objects over 
both fields. Except for one bright z 850 ^b — 24.2 object which met our z-dropout selection 
cut, i-dropouts in our sample range in magnitude from 25.6 down to our completeness limit 
Z850,ab ~ 27.3. For reference, an average-color (UV power-law index ft = —1.3) object 
from Steidel et al. (1999)'s z ~ 3 sample would have ^850,ab ~ 26.1, suggesting a population 
of objects with typical luminosities ranging from ~0.3 to ~1.5 L*. 17 Binning these objects 
into 0.5 magnitude intervals, we illustrate in Figure 6 how the surface density of 2-dropouts 
varies as a function of magnitude. Our typical z-dropout has a half-light radius of 0.15" or 
0.9 kpc, though we find them at all sizes ranging from the limit of the PSF (0.09") to 0.29", 
above which our sample starts to become incomplete. We list all objects which lie in our 
i-dropout sample in Table 1, providing positions, magnitudes, colors, half-light radii, and 
the SExtractor stellarity parameter. Only 13% (~ 1 — 2 objects) of the (z 775 — z$5o) ab > 1-5 
objects without IR magnitudes are likely to be low-redshift contaminants (§2.3). 



4. Predictions 

Before getting into a detailed discussion of the luminosity density at high redshift, 
perhaps the simplest way to begin interpreting what we see at z ~ 6 is to compare it with 
the z ~ 5 ^-dropout sample we previously selected from the HDF North and South (BBI). 
In that work, we used our cloning machinery (Bouwens, Broadhurst & Silk 1998a,b; BBI) 
to project the z ~ 3 £7-dropout population to z ~ 4 — 5 for comparison with our z~4S 
and z ~ 5 V-dropout samples. We found an overall drop in the rest-frame continuum UV 
luminosity density (46% decrease to z ~ 5) as well as a modest decrease in the physical size 
of objects at higher redshifts relative to the z ~ 3 population. 

We would now like to move this comparison out to higher redshift, using our HDF 
^-dropout sample to make an estimate for the surface density of z-dropouts on the sky. 
Unfortunately, the HDF North ^-dropout sample we derived in that work ((V 60 6 — I&ia)ab > 
1.5, (Veoe - Isu)ab > 3.8, (I 8 u - H im ) AB - 1.54, I 8 u,ab > 24, and I 8 u,ab < 27.6) only 



17 Our Z850.AB = 24.2 object (1252-5224-4599) would therefore be a rather bright and suspiciously rare 
~ 6L* object. Even if this object is a low-redshift contaminant (and we have no reason to believe that it 
is), it does not appear to represent a very large source of contamination, given the lack of similarly bright 
objects to zab ~ 25.6 and relative homogeneity of objects faintward of that. 



-16- 




Fig. 6. — Comparison of the number counts for the z-dropouts (i — z > 1.5) observed in our 
two fields (histogram), corrected by 13% for possible contamination, with the no-evolution 
expectations based upon our GOODS z ~ 5 l^-dropout sample (shaded regions). We also 
include the predictions for two different imaging depths: three orbits in z 775 /five orbits in 
2 8 5o (solid line) and six orbits in 2 775 /ten orbits in 2; 850 (dotted line). This illustrates how 
important the issue of completeness is at these depths. (Note that the shaded region above 
assume that 65% of our selection area was at the shallower of these two depths and 35% of 
the area was at the deeper.) For reference, the thin and thick dashed lines show the z ~ 3 
Steidel et al. (1999) LF placed at z ~ 6 (m z ,* = 26.1, a = -1.6, <f) = 0.0025 Mpc" 3 ) with no 
change in normalization (thin) and a normalization that is 4x lower (thick), respectively. (An 
i-dropout selection volume of 800 Mpc 3 arcmin -2 is assumed for both.) The observations, 
while slightly in excess of the z ~ 5 no-evolution predictions (23±25%), are consistent with 
no significant evolution over the redshift interval z ~ 5 — 6. 
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contains 1 object brighter than Iab = 26.5 and hence is not extremely useful in this regard. 
(The fainter ^-dropouts would not be detectable at z ~ 6 to the current depths.) We 
therefore resort to use of a much larger-area ^-dropout sample from the CDF South GOODS 
fields (see Bouwens et al. 2003, in preparation). The virtue of this sample is both its size (130 
objects) and its ACS (PSF FWHM ~0.08") resolution. Our selection criteria for this sample 
were {V m Q-i 77b ) A B > 1-7, (V m6 -i 775 ) A B > 1.1875(z' 7 75-^85o)^s+l-225, (i 775 -z 8 5o)AB < 1-2, 
z 850,ab < 27.2, where the SExtractor stellarity was less than 0.85 in the z 850 image (non- 
stellar with high confidence). We selected this sample from essentially the entire area of 
the CDF South, i.e., 150 arcmin 2 (~30x the area of the HDF), so it should be moderately 
representative of the universe at z ~ 5. Compared to our HDF sample, this ^-dropout 
sample has a similar distribution of redshifts and luminosities. The methodology for the 
generation of this sample remains that of BBI. 

With our cloning machinery, we project objects from our Ve;o6-dropout sample to higher 
redshift using the product of the volume density, l/V max and the cosmological volume. We 
explicitly include pixel-by-pixel k-corrections, cosmic surface brightness dimming, and PSF 
variations in calculating the appearance of these objects to z ~ 5.5 and beyond. Assuming 
no-evolution in the properties of these samples, we ran Monte-Carlo simulations to predict the 
number of i-dropouts that would be observed in our RDCS1252-2927 field and HDF North 
fields. We explored searches at two different depths: three orbits in i 775 , five orbits in z 850 ver- 
sus six orbits in i 775 , ten orbits in z$5o- We then weighted these predictions by the fractional 
area observed to these two different depths (65% at the shallower depth, 35% at the deeper). 
The result is that we predict finding 12.1±1.0 and 4.7±0.4 z-dropouts in our RDCS1252-2927 
and HDF North fields, respectively. The quoted uncertainties reflect the finite size of the 
input V606-dropout sample (see Bouwens et al. 1998a,b). Adding sample variance (simple 
Poissonian statistics), this works out to 12.1±3.6 2-dropouts in our RDCS1252-2927 field 
versus the 20 observed (18.7 after the 13% correction for the low-redshift contamination) 
and 4.7±2.2 z-dropouts in our HDF North field versus the 2 observed (1.9 after correction 
for contamination), indicating a slight reduction in the numbers from z ~ 6 to z ~ 5, but 
more realistically consistent with no evolution (see Table 2). We combine the two fields to 
derive z-dropout number counts and again compare with the no-evolution predictions (Fig- 
ure 6), finding 23±25% more z-dropouts than are predicted using our no-evolution model. 
(Assuming that the luminosity density at z ~ 6 is proportional to the light in the integrated 
counts f(10~ 0Am )^dm, we infer that the luminosity density is 51±29% higher at z ~ 6 
than it is at z ~ 5. Removing the bright zab ~ 24.2 object lowers the luminosity increase 
to 18±23%.) We include on this figure the predictions for the two depths described above 
(solid and dotted lines), showing the effect the assumed depth can have on the predicted 
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numbers. 18 Taking the apparent completeness into account (see also Figure 8 and the discus- 
sion in §5), we quote an approximate surface density of ~0.5±0.2 i-dropouts arcmin -2 down 
to zab ~ 26.5. For reference, we also include a comparison of the predicted and observed 
redshift distributions (Figure 7). 



5. Estimated Luminosity Density 

In the previous section, we compared the z-dropouts we observe with a no-evolution 
projection of a wide-area z ~ 5 V<306-dropout sample. Not only do these simulations give us 
some gauge of the evolution across this redshift interval, but they provide an approximate 
estimate for the i-dropout selection function (assuming no-evolution in the size, shape, or 
color distribution). It is now relatively straightforward to make three different estimates of 
the z ~ 6 luminosity density. 

Our first and most direct estimate comes directly from the comparisons presented in §4, 
where marginal evolution is observed from z ~ 6 to z ~ 5. Linking this with the result from 
BBI (46% drop from z ~ 3 to z ~ 5) yields a 20 ± 29% decrease in the luminosity density 
from z ~ 3 to z ~ 6 (a 37% decrease ignoring the bright zab ~ 24.2 object). This relative 
decrease can then, in turn, be expressed as an absolute luminosity density integrating the 
Steidel et al. (2000) LF down to zab ~ 27 or ~ 0.5L*, where our z-dropouts counts are 
clearly becoming incomplete. The result is 9.8 ± 3.6 x 10 25 ergss _1 Hz -1 Mpc -3 . 

Our second estimate closely follows the more standard approach pioneered by Steidel 
et al. (1999) in deriving the luminosity density at z ~ 3. With this approach, one derives a 
C/y-continuum luminosity function <fi{M) as follows: 

4>{m) — N{m) 



Veff(m) 



'eff 

where M is the absolute magnitude at 1600 A corresponding to some z^q magnitude m 
assuming a fixed redshift of 5.9 (the average redshift for an i-dropout, see Figure 7). The 
effective volume then is calculated as a function of magnitude V e ff(m) = J z p(m, z)^-dz, 
where p(m, z) is the probability that an object of magnitude m and redshift z falls in our i- 
dropout sample and ^ is the cosmological volume at redshift z. The factor p(m, z) contains 



18 We observe a similarly strong dependence in the data, finding 65% more i-dropouts arcmin -2 in the 
deeper overlap regions of our RDCSi252-2927 field than in regions at just half that depth (0.4™ deeper to 
the same S/N). This illustrates how important a consideration incompleteness can be in the magnitude range 
we are considering. 
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Fig. 7. — Comparison of the estimated redshift distribution for the observed z-dropouts 
versus that predicted based upon our GOODS z ~ 5 V-dropout sample (shaded regions). 
This shows that the bulk of the sample is expected to lie between z ~ 5.5 and z ~ 6.2. 
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a whole range of different selection effects which affect the inclusion of an object in our sample 
from the intrinsic distribution of colors to photometric scatter to the effect of the inherent 
surface brightnesses on the completeness of the sample. Assuming a similar distribution 
of surface brightnesses, shapes, and colors to that seen at z ~ 5, we can use the Monte- 
Carlo cloning simulations presented in §4 to determine this function. (Due to the lack of 
V(j06 dropouts to probe the selection function p(m, z) at bright magnitudes, we require that 
p(m, z) be a strictly decreasing function of magnitude.) We present our result in Figure 8. 

For this approach to be effective, the absolute magnitude M has to be a very tight 
function of apparent magnitude m. This occurs when the selection function p(m, z) is a 
very narrow function of redshift. Since this is not the case, we rewrite the above expressions 
as J 0(M(m, z))p(m, z)^dz = N(m) where the absolute magnitude M is a function of 
the apparent magnitude m and the redshift z. 19 Approximating 4>(M) and N(m) as a 
series of step functions, we are able to invert Eq. (3) to solve for (f)(M) (see Figure 9 
for a comparison of the derived LF with the Steidel et al. f999 z ~ 3 determination). 
Integrating this LF down to zab ~ 27 (0.5L*), we find 7.2 ± 2.5 x 10 25 ergss _1 Hz -1 Mpc~ 3 
(5.9 ± 1.8 x 10 25 ergss" 1 Hz -1 Mpc -3 ignoring the bright zab ~ 24.2 object). This represents 
a 39 ± 21% drop relative to what Steidel et al. (1999) report at z ~ 3 to a similar limiting 
luminosity. Note that for z ~ 6 L*-type objects, the effective survey volume is approximately 
3 x 10 4 Mpc. 

In the first two approaches presented, there is the implicit assumption that the se- 
lectability of z ~ 6 objects is similar to that found at z ~ 5, both in their colors and their 
surface brightnesses, and to a large extent, this is probably true. The distribution of z^q — J 
colors (while perhaps a little bluer) are not that different from the low-redshift expectations 
(Steidel et al. 1999; BBI) and similarly for the distribution of surface brightnesses (as similar- 
ities between the predicted and observed number counts in Figure 6 effectively illustrates). 20 
This being said, given the expectation that higher redshift objects are denser and therefore 
of higher surface brightness, it is useful to consider a third approach where no completeness 
correction is made. 



19 In order to convert the apparent magnitudes we measure to absolute magnitudes, it was necessary for 
us to take into consideration possible biases in the measurement of the zsso magnitudes. We can make an 
estimate for this bias using the simulations presented in §4. Comparing the magnitudes we recover to those 
expected based upon an extrapolation of the original ^-dropout photometry to z ~ 5.5, we find an average 
~0.2 m faintward offset in the magnitudes. 

20 Dctermining both distributions observationally would require the sort of deep optical and infrared images 
promised by the HST Ultra Deep Field. Lacking such, one is forced to assume a certain similarity to lower 
redshift dropout populations (as we have done here). 
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Fig. 8. — The probability p(m, z) of some object of z 850j ab magnitude m and redshift z 
being included in our i-dropout sample. This function was computed from a no-evolution 
projection of our wide area V606-dropout sample to z ~ 6 (§4). 
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Fig. 9.— The z-dropout LF derived from our RDCS1252-2927 and HDF North fields (solid 
dots) using a generalized version of Steidel et al. (1999)'s V e ff(m) technique. The z ~ 3 
LF of Steidel et al. (1999) is superimposed as a solid line (Mi6oo,ab — A^i7oo,ab = 0.14 was 
assumed.) We also include this same LF with a 39% lower normalization (dotted line) to 
match the observed evolution from z ~ 6. The error bars represent la uncertainties. 
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This enables us to put a lower bound on the z ~ 6 space density in the event that these 
changes have a substantial effect on the completeness of the 2-dropout population. As with 
the previous two approaches, we make use of the simulations from §4. The difference is that 
we only consider objects which actually make it into our z-dropout sample when computing 
V e ff(m), not every object from our simulations. The effective volume, V e ff(m), is hence 
computed as an ensemble average over all inputs observed at magnitude m and selected as 
i-dropouts. Putting in our observed surface densities for these same magnitude intervals (Fig- 
ure 6), we derive a UV-continuum luminosity density of 2.7 ±0.6 x 10 25 ergss -1 Hz -1 Mpc -3 . 
We emphasize that relative to our previous estimates, this final estimate represents a strict 
lower limit on the luminosity density. 

We converted these UV luminosity densities to star formation rate densities using the 
relation 

SFR 

L uv = const x -ergs s -1 Hz -1 (2) 

Moyr" 1 

where const = 8.0 x 10 27 at 1500 A for a Salpeter IMF (Madau et al. 1998). The result is 
0.0123±0.0045 M yr" 1 Mpc" 3 , 0.0090±0.0031 M Q yr" 1 Mpc" 3 , and 0.0034±0.0008 M yr -1 Mpc~ 3 , 
respectively, for the three approaches just presented. Assuming a Schechter luminosity 
function with faint end slope a = —1.6 and extrapolating this to the faint end limit 
yields a SFR density which is ~ 4x larger. This works out to an integrated star for- 
mation rate of ~ 0.049 ± 0.018 M yr -1 Mpc~ 3 , ~ 0.036 ± 0.012 Moyr^Mpc -3 , and ~ 
0.014 ± 0.003 M yr -1 Mpc -3 , respectively, for these three approaches. To put these esti- 
mates in context, we make a comparison with several previous determinations (Steidel et al. 
1999; Madau et al. 1998; Thompson et al. 2001; Lilly et al. 1996; Stanway et al. 2003; BBI) 
in Figure 10, truncating the observationally-derived LFs at similar faint-end luminosities. 

6. Discussion 

In this work, the luminosity density at z ~ 6 is evaluated using three different proce- 
dures. The first is strictly differential in nature. V^-dropouts from the wide-area GOODS 
survey are projected to z ~ 6 for comparison with the observed i-dropouts. Adopting previ- 
ous results to z ~ 5 (BBI; Thompson et al. 2001) then implies a ~40% drop in the luminosity 
density from z ~ 3 to z ~ 6. The second and third procedures, by contrast, are more direct, 
relying on a derived selection function to convert the observed z-dropouts into a luminosity 
density at z ~ 6. Relative to the luminosity density reported at z ~ 3 by Steidel et al. (1999), 
the resulting drop is a factor of ~ 2 and ~ 5, depending upon the two different assumptions 
these two procedures make about the completeness levels (and therefore the presence of low 
surface brightness objects at high redshift). Similarities between our first two estimates of 
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Fig. 10. — A history of the star formation rate density assuming no extinction correction, 
integrated down to 0.5L*. We include three determinations from this work (the large solid 
red circles), the lower point assuming no incompleteness correction and providing a reliable 
lower limit, the slightly larger middle point based on a generalization of Steidel et al. (1999) 's 
V e ff(m) formalism, and the upper point based on the differential evolution measured from 
z ~ 6 to z ~ 5 (and linked to z ~ 3 using the results quoted in BBI) (see §5 for details). 
The middle point gives our preferred estimate. A Salpeter (1955) IMF is used to convert 
the luminosity density into a star formation rate (see, for example, Madau et al. 1998). Our 
different estimates provide a nice illustration of the real uncertainties in the star formation 
rate density at z ~ 6, arising from strong redshift-dependent selection effects. This topic 
will be explored much more extensively in Bouwens et al. (2003). Comparison is made with 
the previous high redshift determinations of Lilly et al. (1996) (open squares), Madau et al. 
(1998) (open circles), Steidel et al. (1999) (crosses), Thompson et al. (2001) (open triangles), 
Stanway et al. (2003) (open pentagons), and BBI (solid red triangles). The top horizontal 



-25 - 



the luminosity density point to a general consistency between those methodologies. 

Our three estimates have different strengths and weaknesses. Our first estimate, for 
example, depends on the relatively small area BBI result (HDF-N and HDF-S) and therefore 
could be quite sensitive to cosmic variance. Our second estimate, by contrast, relies upon 
the much larger area z ~ 5 V-dropout population (CDFS) to derive the z ~ 6 selection 
function. Unfortunately, this same z ~ 5 sample could be subject to similar selection effects 
and therefore missing light, possibly biasing the numbers low. Our first two estimates provide 
accurate estimates of the z ~ 6 luminosity density in lieu of large amounts of size or color 
evolution. Our final estimate, on the other hand, is probably much too low except in cases 
where there has been dramatic amounts of evolution in the sizes of galaxies. It therefore 
serves as a useful lower limit. On balance, we prefer our estimate based on a generalized 
version of the Steidel et al. (1999) V e ff(m) formalism, given the effect of cosmic variance on 
our first estimate and the extreme assumptions present in the third (e.g., given the small 
increase in the universal scale factor from z ~ 6 to z ~ 5, the change in surface brightnesses 
should be at most modest). We will use it throughout the remainder of the discussion. 21 

Before getting into a comparison against previous work, we should take a more detailed 
look at the validity of the sample itself, particularly, with regard to our use of a cluster field. 
Two concerns immediately arise. The first regards the effect of lensing by the cluster on 
z ~ 6 sources. We can make a simple estimate of this effect using the Broadhurst, Taylor, 
& Peacock (1995) equation for the magnification bias, M 2,5/i_1 , where \i is the slope of the 
number counts ^l°fLffIrg) _ p or a faint end slope m = 0.24 (a simple consequence of a z ~ 6 LF 
where a = —1.6) and an average magnification of ~1.03, this works out to an expected 1% 
drop in the surface density of i-dropouts relative to the unlensed case. The second concern 
involves the possibility that z ~ 1.2 cluster ellipticals might get mixed up in our i-dropout 
sample. After all, we decided to observe RDCS1252-2927 with the 2775 and zgso filters because 
they straddled the 4000 A break. Fortunately, the mean i — z color for bright ellipticals in 
this cluster is ~ 0.9 (one can see the cluster as a small overdensity in Figure 1 at i — z=0.9 
and z — J =1.1) and this colour decreases quite markedly towards faint luminosities (Postman 
et al. 2002; Blakeslee et al. 2003b). There is no indication that lower luminosity z ~ 1.2 
cluster members are being scattered into the z-dropout sample. 

It is interesting to compare the present results with other early work on the surface 
density of i-dropouts. Stanway et al. (2003), for example, report finding 8 objects over 146 



21 Ideally speaking, we would estimate the z ~ 6 luminosity density using a completely differential proce- 
dure, scaling the sizes and colors to match the evolution observed. Such an analysis is being performed in a 
future paper on the dropouts in the GOODS fields where the samples are substantially larger (Bouwens et 
al. 2003). 
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arcmin 2 , e.g., 0.05 z-dropouts arcmin -2 to zab ~ 25.6. We only find one object (~ 0.02 i- 
dropouts arcmin -2 ) to that same depth, a result which may not be surprising given the 
numbers and relative areas. In a deeper but smaller area search, Yan et al. (2003) report 
finding 2.3 objects arcmin -2 (27 objects) to zab ~ 28.0, subtracting at most 4 objects from 
this estimate due to contamination by stars. Since the median magnitude of their sources is 
z s50,ab ~ 27.4 and their brightest reported source has z 850t AB ~ 26.8, this works out to ~ 1.2 
objects arcmin -2 to our quoted magnitude limit, somewhat larger than the 0.5±0.1 objects 
arcmin -2 we find. We are somewhat uncertain about the reliability of their identifying %- 
dropouts to z S 5o,ab ~ 28 (7 a) in a 5-orbit (9540 s) zgso-band exposure, given the difficulty 
we had in identifying i-dropouts to a much brighter limit of zab ~ 27.3 (6a) in the deeper 
(~10-20 orbit) overlapping regions. For typical z-dropout sizes (~ 0.15"), our analysis would 
tend to suggest a S/N closer to 3 at their stated magnitude limit. 

Compared to Stanway et al. (1995), who estimate a 25% contamination fraction due to 
ellipticals or EROs in an (i — z) > 1.5 sample, we estimate a smaller contamination fraction, 
13%, based upon the measured z — J colors for our sample. Though hardly significant given 
the number of objects involved, such a difference clearly follows the trend toward lower early- 
type fractions at fainter magnitudes. We noted similarly small contamination fractions in 
our analysis of the ^-dropouts from the HDF North (BBI). 

Down to our magnitude limit Zab ~ 27.3, we find a star formation rate density at 
z ~ 6 of 0.0090 ± 0.0031 M yr -1 Mpc -3 , ~2x lower than similar estimates at z ~ 3 from 
Steidel et al. (2003) and Madau et al. (1998). This estimate is nearly 14 x Stanway et al. 
(2003)'s quoted result (6.7 ± 2.7 x 10 -4 M yr -1 Mpc -3 ). While a portion of this difference 
can be attributed to the different depths of our searches (the present census going ~1.4 
magnitudes further down the luminosity function), there are significant differences in the 
effective volumes we assume, the present study correcting for the significant incompleteness 
at z ~ 6 resulting from surface brightness dimming while Stanway et al. (2003) do not make 
such a correction. For both this reason and our greater overall depth, the present study 
represents a clear improvement over the Stanway et al. (2003) estimate. From z ~ 5 to 
z ~ 6, we find minimal evidence for dramatic evolution, consistent with the small amount of 
cosmic time available across this interval (~0.2 Gyr from z ~ 6 to z ~ 5 and ~1 Gyr from 
z ~ 6 to z ~ 3). In general, the trends we find are consistent with the gradual decline in 
star formation density previously reported from z ~ 3 to z ~ 5 (Madau et al. 1996; Steidel 
et al. 1999; Thompson et al. 2001; Lehnert et al. 2003; Deltorn et al. 2003; BBI) or as found 
in recent work on Lyman-alpha emitters, where a deficit is claimed at z ~ 5.7 relative to the 
z ~ 3 population (Maier et al. 2003). 

Since recent work from absorption line studies on z > 5.8 QSOs indicates that the 
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universe may have been reionized near a redshift of z ~ 6 (Fan et al. 2001; Becker et al. 
2001; Fan et al. 2002), it is relevant to compare the star formation rate we determine here 
with that needed to reionize the universe at z ~ 6. For the latter, we use the following relation 
from Madau, Haardt, & Rees (1999), correcting it for the baryon density [Vt^h 2 = 0.0224) 
derived from the recent WMAP results (Bennett et al. 2003) and shifting it to z ~ 6: 

p\ « (0.052 M yr" 1 Mpc -3 ) (^) C 30 (±±*\ ' . (3) 



esc • 



where p„ is the star formation rate density, C30 is the H 1 concentration factor (1/30) ^p|^ i 

{PH 1 ) 2 ' anc ^ ^ esc * s ^ ne f rac ti° n °f ionizing radiation escaping into the intergalactic medium 
Given the large observational and theoretical uncertainties in the exact values of C30 and /, 
the star formation rate densities inferred here at z ~ 6 (~O.OO9OM yr -1 Mpc -3 observed, 
~0.036M Q yr -1 Mpc -3 extrapolating the LF to the faint limit) are well within the range of 
that needed to reionize the universe at z ~ 6 (~30% lower than the fiducial value given 
in Eq. (3)). Whether or not the objects we observe at z ~ 6 are sufficient to reionize the 
universe or this ionizing radiation is provided by a completely different population of objects 
(e.g., Madau 1998), it seems clear that this ionizing radiation does not come from z ~ 5 — 6 
AGNs (Haiman, Madau, & Loeb 1998). 



7. Summary 

We use deep ACS WFC i+z observations of RDCS 1252-2927 and the HDF North to 
search for z ~ 6 candidates by looking for a strong Lyman break across the i + z passbands 
in fields with deep IR imaging. We augment this with deep infrared imaging to derive z — J 
colors for the z ~ 6 candidates to help distinguish them from lower redshift interlopers 
where possible. We compare our findings with a no-evolution projection of a wide-area Vqqq- 
dropout sample to z ~ 6 and then use these simulations to make three different estimates of 
the rest-frame continuum C/y-luminosity density at z ~ 6. 

• To zab ~ 27.3, we find 21 2-dropouts (i — z > 1.5; 6cr detections) over 36 arcmin 2 
in our RDCS1252-2927 field (one object is lensed) and 2 z-dropouts over 10 arcmin 2 
in our HDF North field. This is equivalent to ~ 0.5 ±0.1 object arcmin -2 down to 
our magnitude limit zab ~ 27.3, or ~ 0.5 ± 0.2 object arcmin -2 down to zab ~ 26.5 
corrected for completeness. 

• Compared to an average-color (UV power-law index (3 = —1.3) [/-dropout in the 
Steidel et al. (1999) z ~ 3 sample, the z-dropouts we find range in luminosity from 
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~0.3 L* to ~1.5 L* (with the exception of one very bright z$5o,ab = 24.2 (6L*) object 
which also meets our selection criterion). Our typical z-dropout has a half-light radius 
ranging from 0.09" to 0.29", or 0.5 kpc to 1.7 kpc, for a Q M = 0.3, Q A = 0.7, h = 0.7 
cosmology. 

• Using our deep infrared data to constrain the spectral slopes redward of the break, we 
find that all 12 of our z-dropout candidates (i — z > 1.5), for which we have IR data, 
have blue (< 0.8) z — J colors. Using the Kaplan-Maier estimator on all extended 
i — z > 1.3 sources, we estimate that only 13% of such i — z > 1.5 sources are low- 
redshift interlopers. Assuming spectroscopic confirmation of the two-color dropout 
technique to z ~ 6 (e.g., Fan et al. 2001; Bunker et al. 2003), this demonstrates that 
an i — z selection can be an effective way of uncovering a population of high redshift 
objects (see also Stanway et al. 2003). 

• Over the 21 arcmin 2 where we had infrared coverage, we identified only one red (i — z> 
1.3) point-like object with z — J colors consistent with being az~ 5.5 — 6 AGN. Since 
in total we identify 7 point-like objects with zab > 25 and (2775 — z$5o)ab > 1-3 and 
only one out of four objects (25%) with IR coverage has z — J colors less than 0.8, this 
works out to an estimated surface density of ~0.04±0.03 z ~ 5.5 — 6.3 AGNs arcmin -2 
down to our magnitude limit of z^b ~ 27.3. 

• Comparing the number of z-dropouts with a no-evolution projection of our z ~ 5 
V606-dropout sample from CDF South GOODS, we estimate the evolution in rest- 
frame continuum UV luminosity density from z ~ 6 to z ~ 5. We find 23±25% 
more i-dropouts (~21) than are predicted (17) once a correction for the contamination 
rate is made, consistent with no strong evolution over this redshift interval. Redoing 
this increase in terms of integrated luminosity, this increase amounts to 51 ±29% (or 
18 ± 23% removing the one very bright zab ~ 24.2 object). Adopting previous results 
to z ~ 5 (BBI; Thompson et al. 2001), this works out to a mere 20±29% drop in the 
luminosity density from z ~ 3 to z ~ 6. 

• Using simulations based on a set of CDF South V606-dropouts, we estimate the selection 
function p(m, z) for our z-dropout sample. Then, via a generalized version of Steidel 
et al. (1999) 's V e ff(m) formalism, we calculate the [/U-luminosity function for this 
sample, and integrate it down to zab ~ 27.0 (~ 0.5L*). The rest-frame [/U-luminosity 
density we derive (7.2 ± 2.5 x 10 25 ergs s" 1 Hz -1 Mpc -3 ) is 39 ± 21% lower than Steidel 
et al. (1999) found to a similar limiting luminosity, consistent with the above estimate. 
This is our preferred estimate. 



-29- 



• The previous two approaches assume no large change in the selectability (color or sur- 
face brightness distribution) of dropouts from z ~ 5 to z ~ 6. Given the expectations 
(observational and theoretical) that dropouts may have higher surface brightnesses at 
z ~ 6 than z ~ 5, it is useful to make a third estimate for the luminosity density 
from the simulations, but this time assuming no incompleteness. Running through the 
numbers, we find 2.7 ±0.6 x 10 25 ergss -1 Hz -1 Mpc~ 3 for the observed i-dropouts, ~5x 
lower than the value reported by Steidel et al. (1999) to a similar luminosity. This 
third approach provides a reliable lower limit to the luminosity density. 

• Converting the luminosities densities we infer into star formation rate densities using 
standard assumptions (e.g., Madau et al. 1998), we find an integrated star forma- 
tion rate density of 0.0123 ± 0.0045 M yr -1 Mpc -3 , 0.0090 ± 0.0031 M yr -1 Mpc -3 , 
and 0.0034 ± 0.0008 M Q yr -1 Mpc -3 , respectively, for the three approaches just pre- 
sented (or ~ 0.049 ± 0.018 M yr -1 Mpc -3 , ~ 0.036 ± 0.012 M yr -1 Mpc -3 , and 
~ 0.014 ± 0.003 M yr -1 Mpc -3 extrapolating the observations to low luminosities 
using a Schechter function with faint-end slope a = —1.6). 

• Our preferred estimate for the rest-frame continuum UV luminosity density and star 
formation rate density at z ~ 6 are 7.2 ± 2.5 x 10 25 ergss -1 Hz -1 Mpc -3 and 0.0090 ± 
0.0031 M yr -1 Mpc" 3 (~ 0.036 ± 0.012 M Q yr -1 Mpc" 3 extrapolating the luminosity 
function to the faint limit), respectively. This represents a 39±21% drop from z ~ 3 
to z ~ 6. 

• The z ~ 6 rest-frame continuum ?7U-luminosity densities we infer are well within the 
expected range needed for reionization, for canonical assumptions about the H i clump- 
ing factor and the fraction of UV radiation escaping into the intergalactic medium. 

The combination of deep ACS % + z and ground-based IR imaging have been shown to 
be a very effective means of isolating high redshift objects and studying their properties. We 
will be following up this analysis with an investigation of the shallower, but larger area, data 
in the CDF South and HDF North from the GOODS program (Bouwens et al. 2003). 

We extend a special thanks to Mark Dickinson for providing us with his fully reduced 
NICMOS images of the HDF North. We would also like to acknowledge T. Allen, K. An- 
derson, S. Barkhouser, S. Busching, A. Framarini, and W.J. McCann for their invaluable 
contributions to the ACS Investigation Definition Team (IDT). This research has made use 
of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under contract with the National Aeronautics 
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and Space Administration. ACS was developed under NASA contract NAS 5-32865. RJB, 
GDI, and the ACS IDT acknowledge the support of NASA grant NAG5-7697. 
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This preprint was prepared with the AAS IATgX macros v5.0. 
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A. Noise 

The drizzling and resampling procedure we employ (and generally employed with HST 
data sets) to produce our fully coaligned optical and infrared data set naturally introduces a 
certain correlation into the noise of this data set. To estimate the true amplitude of the noise- 
distinct from the single pixel RMS-as well as the effective noise kernel, RMS fluctuations 
within apertures of increasing sizes are measured, after masking out > 5a detections, up to 
a scale where there is no obvious additional power. A noise model (amplitude and kernel) 
is then constructed which provides a plausible fit to the observations. This model is used 
throughout to estimate errors in the photometry. 



B. z ~ 6 AGN 

We found one (i 775 — Zs5o)ab ~ 1.3, (zsso — J)ab < —0.1 stellar object consistent 
with a z ~ 5.5 — 6 AGN-identification over the area where we had infrared coverage (~21 
arcmin 2 ) (Figure 11). Since in total we identify 7 point-like objects with z^b > 25 and 
(i 775 — z 850 )ab > 1-3 and only one out of the four objects with IR coverage (25%) had 
(^850 — J)ab colors less than 0.8, this works out to an estimated surface density of ~0.04±0.03 
objects arcmin -2 down to our magnitude limit of z 850) ab ~ 27.3. 
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Fig. 11. — ^775^850 JiiT images of a z ~ 6 AGN candidate found in our RDCS1252-2927 field. 
The above object is the only such candidate over the 21 arcmin 2 where we have both optical 
and infrared imaging. The position in % — z, z — J color-color space is included along with a 
plausible fit to an SED. 
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Table 1. z ~ 6 Sample. Objects from RDCS1252-2927 and the HDF North which satisfy 

our (z 775 — z$5o)ab > 1-5 z-dropout criterion. 



Object ID Right Ascension Declination Z850 a i — z° z — J c V — z I — z S/G d rhi(") 



1252-5224-4599 6 


12 


:52 


:56 


.8880 


-29 


:25 


:55 


.503 


24 


,2±0. 


1 


1 


.5 





1 









0. 


29 


1252-2134-1498 


12 


:52 


:45 


.3816 


-29: 


:28: 


27 


.105 


25 


,6±0. 


2 


>2 


,1 





,4 






0.01 


0. 


18 


1252-2585-3351 


12 


:52 


:52 


.2832 


-29: 


:28: 


:04 


.743 


25 


,7±0, 


1 


2 


.0 


0, 


1 






o 


0. 


20 


1252-6031-966 


12 




:43 


.3793 


-29 


:25: 


17 


.474 


25 


,7±0, 


1 


2 


.0 










0.30 


0. 


11 


1252-562-2836 


12 


:52 


:50 


.2971 


-29 


:29: 


43 


.226 


25 


,8±0. 


1 


2 


.0 










0.06 


0. 


.14 


1252-6798-302 


12 


:52 


:40 


.8435 


-29: 


:24: 


39 


.033 


25 


,8±0, 


2 


1 


.7 










0.04 


0. 


16 


1252-4720-6466 


12 


.53 


:03 


.8342 


-29: 


.26: 


20 


.835 


25 


,8±0, 


1 


1 


.5 




- 


- 


- 


0.02 


0. 


12 


1252-5058-5920 


12 


:53 


:01 


.7448 


-29 


:26: 


03 


.873 


25 


9±0, 


2 


>1 


.8 


-0. 


4 









0. 


19 


1252-277-3385 


12 


:52 


:52 


.3987 


-29: 


:29: 


57 


.525 


26 


,0±0, 


2 


>1 


.6 













0. 


19 


1252-5377-2621 


12 


:52 


:49 


.5122 


-29: 


:25: 


:47 


.721 


26 


0±0, 


2 


>1 


.6 


<-0 


.1 









0. 


22 


1252-7065-6877 


12 


:53 


:05 


.4237 


-29: 


:24: 


26 


.222 


26 


,1±0, 


2 


1 


.6 










0.01 


0. 


.14 


1252-3544-3542* 


12 


:52 


:53 


.0237 


-29 


:27 


16 


.801 


26 


2±0, 


2 


1 


.6 





.1 









0. 


23 


1252-2439-3364 


12 


:52 


:52 


.3340 


-29: 


:28: 


12 


.030 


26 


,4±0, 


2 


>1 


.7 


<-0 


2 









0. 


17 


1252-7005-1697 


12 


:52 


:46 


.1876 


-29: 


:24 


:28 


.846 


26 


5±0, 


2 


>1 


.5 













0. 


13 


1252-3729-4565 


12 


:52 


:56 


.7460 


-29: 


:27 


:07 


.631 


26 


,7±0, 


2 


>1 


.8 


0, 


.4 









0. 


12 


1252-7313-6944 


12 


:53 


:05 


.6814 


-29: 


:24 


13 


.836 


26 


,7±0, 


2 


>1 


.6 










0.09 


0. 


09 


1252-1199-3650 


12 


:52 


:53 


.4200 


-29: 


:29 


11 


.443 


26 


,7±0, 


2 


>1 


.7 





.1 









0. 


12 


1252-5399-4314 


12 


:52 


:55 


.7957 


-29: 


:25 


:46 


.724 


27 


,0±0, 


2 


>1 


.6 


<0 


.0 









0. 


11 


1252-3497-809 


12 


:52 


:42 


.7537 


-29 


:27 


:18 


.894 


27 


0±0, 


2 


>1 


.6 













0. 


11 


1252-4375-1877 


12 


:52 


:4G 


.8536 


-29 


:2G 


:37 


.733 


27 


2±0, 


2 


1 


.6 


<0 


.7 






0.60 


0. 


09 


1252-3684-528 


12 


:52 


:41 


.6802 


-29: 


:27: 


:09 


.550 


27 


,3±0, 


2 


>1 


.5 










0.16 


0. 


09 


HDFN-2135-3286 


12 


.36 


:49 


.9368 


62 


:13: 


:55 


.671 


26 


5±0, 


2 


>1 


.9 


<-0 




>2.1 


1.1 


0.02 


0. 


14 


HDFN-4965-4355 


12 


:36 


:30 


.7457 


62 


12 


:53 


.371 


27 


,0±0, 


2 


>1 


.5 






>1.3 




0.01 


0. 


.10 



a AB Magnitudes 
b All limits are 2a. 

c Hcre the J band alternatively refers to the ISAAC J band and the NICMOS Jn filter depending upon which field 
the object is found. 

d SExtractor stellarity parameter, for which = extended object and 1 = point source 
c The "1252-" prefix denotes an object from RDCS1252-2927. 
*This object appears to be lensed and therefore is excluded from our sample. 
**This is {z-J U o)ab- 



-36- 



Table 2. The number of z-dropouts found in our samples versus no-evolution predictions. 
Two different la uncertainties are quoted on all predictions, the first based on the finite 
size of our ^-dropout sample (130 objects) (Bouwens et al. 1998a,b) and the second based 

on sample variance (simple Poissonian errors). 



Data set 


Observed 


No evolution Prediction 


RDCS1252-2927 


18.7 a 


12.1±1.0±3.5 


HDF North 


1.9 a 


4.7±0.4±2.2 



a Corrected for the expected 13% contamination rate. 



